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1 Introduction 

Some one-dinicnsional quantum spin chain Hamiltonians and some models of classical statistical mechan- 
ics in two spatial dimensions on a lattice - the vertex models - share a common mathematical structure 
responsible by well-known breakthroughs in our present understanding of these models [1, 2, 3]. If the 
Boltzmann weights underlying the vertex models are obtained from solutions of the Yang-Baxter (YB) 
equation, the commutativity of a collection of transfer matrices depending on a spectral parameter, say 
u, follows, leading to their integrability. Moreover, by taking the logarithmic derivative of the transfer 
matrix, evaluated at a special value of this spectral parameter, one gets an associated one-dimensional 
quantum spin chain Hamiltonian. If fact, it might be said that the most successful approach to construct 
integrable two-dimensional lattice models of statistical mechanics is by solving the Yang-Baxter equation. 
For a given solution of this equation, one can define local Boltzmann weights to find a commuting family 
of transfer matrices. 

The diagonalization of one-dimensional quantum spin chain Hamiltonians started with the Bethe 
Ansatz (BA) [4], which gradually become a powerful method in the analysis of integrable models. There 
are several versions: the coordinate BA [4], the algebraic BA [5], the analytical BA [6], etc. The simplest 
version is the coordinate BA. In this framework one can obtain the eigenfunctions and the spectrum of 
the Hamiltonian from its eigenvalue problem. It is really simple and clear for two-state models like the 
six-vertex model but becomes tricky for models with a higher number of states. 

The development of the quantum inverse scattering method approach to integrability resulted in the 
algebraic BA [5] which is an elegant and important generalization of the coordinate BA. It is based on 
the idea of constructing the eigenvectors of the transfer matrix via the action of " creation " operators 
on a reference state. The creation operators are just entries of the monodromy matrix whose trace is the 
transfer matrix. Then by using the YB equation one can write the so-called fundamental relation which 
entails the generalized " commutation relations " of the transfer matrix with the remaining entries of the 
monodromy matrix ( the " creation " and " annihilation " operators) . The creation operators that form a 
n-body Bethe vector are evaluated at some unknown spots, say {ui} ~ {ui, U2, ..., u„}. Finally, the action 
of the transfer matrix, evaluated at some spectral parameter, say u, will be an eingenvalue problem, if 
a set of equations, known as Bethe equations, not depending on u, determine the set {ui}. Despite its 
elegance and completeness, the actual implementation of the algebraic Bethe Ansatz can become rather 
tricky and laborious, as will be exemplified in this article. It is usual to rely on other methods such as the 
coordinate or the analytical BA to gather important informations about the eingenvalues of the transfer 
matrix. 

The integrability of open spin chains in the framework of the quantum inverse scattering method was 
formulated by Sklyanin relying on previous results of Cherednik [7]. In reference [8], Sklyanin used his 



formalism to solve the open spin-1/2 chain with diagonal boundary terms, and his original formalism 
was successfully extended to others models by Mezincescu and Nepomechie in [9]. The basic idea is 
similar to that of the usual algebraic BA, with the important difference that now one considers the so- 
called double-row monodromy matrix. When considering systems on a finite interval, with independent 
boundary conditions at each end, one has to introduce reflection matrices to describe such boundary 
conditions. Integrable models with boundaries can be constructed out of a pair of reflection i^-matrices 
K (u) which obey the cousin equation of the YB equation - called reflection equation. K^ (u) and K^{u) 
describe the effects of the presence of boundaries at the left and the right ends, respectively. Sklyanin has 
shown that the double-row monodromy matrix obeys a generalization of the fundamental relation also 
named reflection equation. The non-diagonal entries of the double-row monodromy matrix play the role 
of "creation " and "annihilation " operators and the diagonal ones are used to construct the transfer 
matrix, the reflection equation providing the fundamental set of "commutation relations ". Similarly 
to the usual BA the creation operators acting on a reference state form a n-body Bethe vector when 
evaluated at some unknown spots, say {ui} — {ui, U2, ..., Un}, and an eingenvalue problem for the transfer 
will exist if a set of equations, named also Bethe equations, not depending on u, determine the set {ui\. 
In a recent article [10] Guang-Liang Li, et al investigated the 19- vertex Izergin-Korepin (IK) model [11] 
with boundaries using the algebraic BA. They have found a spurious dependence of the Bethe equations on 
the spectral parameter u , such that the set {ui\ could inconsistently depend on the spectral parameter u. 
Even more, this spurious dependence led the authors of [10] to raise the suspicion of the non-uniqueness 
(or inconsistency) of Skyanin's algebraic BA for this model. On the other hand, this model has been 
previously studied through the analytical BA [12, 13] and the coordinate BA [14] where no such spurious 
dependence was found. So, we decided to investigate a related 19-vertex model , the Zamolodchikov- 
Fateev (ZF) model [15], where one surely does not expect such a spurious dependence on u. The approach 
is the same as that of [10], first envisaged by Fan [16], which is the boundary version of the Tarasov work 
[17] for the periodic IK model. As a further check we also performed the calculations for the IK model . 
We do not have found such a spurious dependence for the ZF model, our results agreeing with all previous 
results. For our pleasant surprise we do not have found such a spurious dependence for the IK model too. 
Indeed, our calculations agree with those of Guang-Liang Li, et al and there is no spurious dependence on 
the spectral parameter at all. The point is that this spurious dependence of the Bethe equations on the 
spectral parameter in [10] is only apparent. A careful analysis as we perform em Section 4 bellow shows 
that it does not exist. It can be also checked that, in the notation of ref. [10], that their factor function 
(3{u,Ui) in their Eq. (97) does not depend on u, so that the dependence on u in their Eq (102), for the 
non-quantum group invariant case, is non existent. The Sklyanin algebraic BA is unique and consistent 
for the ZF and the IK vertex models. 

The paper is organized as follows: in Section 2 we define the models to be studied and introduce the 



necessary notation to briefly review the basic concepts of tlie algebraic BA.In section 3, we present our 
detailed calculations common to both models and in Section 4 the eigenspectra and the corresponding 
Bethe equations are explicitly presented for each model. Section 5 is reserved for conclusions and in 
the appendix we derive the fundamental commutation relations for the double-row monodromy matrix 
entries. 

2 The Models 

To determine an integrable vertex model on a lattice it is first necessary that the bulk vertex weights be 
specified by a 7?.- matrix TZ(u), where u is the spectral parameter. It acts on the tensor product V (!S)V for 
a given vector space V and satisfy a special system of functional equations, the Yang-Baxter equation 

ni2{u)ni3{u + v)n23{v) = 7^23(«)7^l3(u + v)tzi2{u), (2.1) 

inV (^V®V, where 7^l2 = 7?. «> 1, 7^23 = 1 ® 7^, etc. 

The TZ matrix is said to be regular if it satisfies the property TZ{0) = P, where P is the permutation 
matrix in F T^: P{\ct) ® |/3)) — |/3) \a) for \a) , |/3) G V. In addition, we will require [9] that TZ{u) 
satisfies the following properties 

regularity : n^O) ^ f{0y/^Pi2, 

unitarity : TZi2{u)n{'2%-u) ^ f{u), 

PT - symmetry : Pl27^l2(w)Pl2 = Kl^'" («), 

crossing — symmetry : 7?.i2(u) = C/i7?.*|(— u — p)f7f "^, (2.2) 

where f{u) — xi{u)xi{—u), xi(u) being defined for each model bellow, ti denotes transposition in the 
space i , p is the crossing parameter and U determines the crossing matrix 

M^U^U^MK (2.3) 

Unitarity and crossing-symmetry together imply the useful relation 

Ml7^*|(-« - p)Mf l7^*^2(« -p) = f{u). (2.4) 

The boundary weights then follow from iiT-matrices which satisfy the boundary versions of the Yang- 
Baxter equation [8, 9], the refiection equation 

7^l2(u - v)K^{u)n^^^^{u + v)K^{v) ^ K^{v)ni2iu + v)K^{u)n[\*^{u - v), (2.5) 

and the dual reflection equation 

TZui-u + v){K+Y^ {u)M^^n['^''{^u -V- 2p)Mi{K+y^v) 



= {K+Y^{v)MiTZi2{-u - V - 2p)M-\K+Y^{u)TZ%*''{-u + v). (2.6) 

In this case there is an isomorphism between K^ and K^ : 

K-{u):^K+{u)^K-{-u- pfM. (2.7) 

Therefore, given a solution of the reflection equation (2.5) one can also find a solution of the dual reflection 
equation (2.6). 

A quantum-integrable system is characterized by the monodromy matrix Tiu) satisfying the funda- 
mental relation 

R{u - v) [T{u) ® T{v)] ^ [T{v) ® T{u)] R{u - v) (2.8) 

where R{u) is given by R{u) = PTZ{u). 

In the framework of the quantum inverse scattering method, the simplest monodromies have become 
known as C operators, the Lax operators, here defined by Caq{u) = TZaq{u), where the subscript a 
represents the auxiliary space, and q represents the quantum space. The monodromy matrix T{u) is 
defined as the matrix product of N Lax operators on all sites of the lattice, 

T{u) = CaN{u)CaN-l{u) ■ ■ ■ Cal{u). (2.9) 

We recall that the main result for integrability is that, if the boundary equations arc satisfied, then 
the Sklyanin's transfer matrix 

t{u) = Txa {K+{u)T{u)K-{u)T-^{-u)) , (2.10) 

forms a commuting collection of operators in the quantum space 

[i(u),i(w)] ==0, Vu,w (2.11) 

The commutativity of t{u) can be proved by using the unitarity and crossing-unitarity relations, the 
reflection equation and the dual reflection equation. In particular, it implies the integrability of an open 
quantum spin chain whose Hamiltonian (with K^ifi) = 1) is given by 

N-l 



1 dK-{u) 

k,k+l + 77 



du 
fc=i 



tY„K+{0)HN.o 
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where the two-site terms are given by 

Hk,k+1 — -rPk.k+iT^k.k+i{u) , (2-13) 

in the standard fashion. 



The three-state vertex models that we wiU consider are the Zamolodchikov-Fateev model and the 
Izergin-Korcpin model. Their 7?,-matrices have a common form 
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satisfying the properties (2.1-2.4). 

In the context of the coordinate BA, these models were solved in [14], where an appropriate parametriza- 
tion of wavefunctions was used in order to recast the coordinate BA for these three-states models in a 
form as simple as the coordinate BA for the two-state models [18]. 

2.1 The Zamolodchikov-Fateev model 

This is the simplest three-state 19-vertex model [15]. In the Bazhanov [19] and Jimbo [20] classification, 
it is the B\ model or the A\ model in the spin-1 representation [21], due to its construction from the 
six-vertex model by the fusion procedure. The 7^-matrix which satisfies the YB equation (2.1) has the 
form (2.14) with 

x\(u) = sinh(M + ?]) sinh(M + 2ry), 0:2(1*) = sinhusinh(u + 77), 

0:3(1*) = sinhMsinh(M — ry), a:4(u) = sinhusinh(M + 77) + sinhry sinh2ry, 

y^(u) — 2:5 (u) = sinh(u + r\) sinh2?7, yek'A ~ x^iu) = sinh?isinh277, 

yji^u) = 2:7 (u) = sinh 77 sinh 277. (2-15) 

This 7?.-niatrix is regular and unitary, with f{u) ~ x\(u)x\{~u\ P- and T-symmetric and crossing- 
symmetric with M =\ and p = r\. The most general diagonal solution for K^ (u) has been obtained in 
Ref. [22] and is given by 



with 



fcii(") 



K-(u,li) = 



(3 sinh u + 2 cosh u 
(3 sinh u ~ 2 cosh u ' 



fcii(w) 



hai^) 



hai^) 



(3 sinh(M + 77) — 2 cosh(u -I- 77) 
(3 sinh(M — 77) + 2 cosh(7i — 77) ' 



(2.16) 



(2.17) 



where /3 is a free parameter. By the automorphism (2.7) the sohition for K^(u) foUows 

K+{u,a)^K-{-u-p,a)= \ 1 , (2.18) 

V k+,{u) J 

with 

_i_ a sinh(u + 77) — 2 cosh(u + 77) , + / \ a sinh u + 2 cosh m 



fciiH = - „,,■ u/.. , ^^ , o u/.. , ^x ' ^33(") = - 



a sinh(M + 77) + 2 cosh(M + 77) ' a sinh(M + 27;) — 2 cosh(7i + 2?7) ' (o to) 

where a is another free parameter. 

The energy spectrum of this model was aheady obtained by Mezincescu at al. [22] through a general- 
ization of the quantum inverse scattering method developed by Sklyanin [8] , the so-called fusion procedure 
[21] and by Fireman et al [14] through a generalization of the coordinate BA. 

For a particular choice of boundary terms, the ZF spin chain has the quantum group symmetry i.e., 
if we choose (3 — 2 coth^_ and a — 2 coth^4-, with ^^ — > cxd, then the open spin chain Hamiltonian has 
C/g(s7i(2))-invariancc [22]. 

The fusion procedure was also used by Yung and Batchclor [13] to solve the ZF vertex-model with 
inhomogeneities. The coordinate and algebraic BA with periodic boundary conditions were presented in 
[23]. 

2.2 The Izergin-Korepin model 

The solution of the YB equation corresponding to Ag in the fundamental representation was found by 
Izergin and Korepin [11]. The 7?,-matrix has the form (2.14) with non-zero entries 

xi{u) = sinh(7i — 577) + sinh77, 2:2(7^) = sinh(M — 377) + sinh377, 

x^{u) = sinh(7i — 7;) + sinh77, a;4(7i) = sinh(u — 377) — sinh577 + sinh377 + sinh?; 

xr,{u) ^ -2e""/2sinh277Cosh(^ - 877), y^{u) = -2c"/2sinh277Cosh(| - 877) 

x^{u) = 2e-"/2+2''sinh277sinh(|), 776(w) = -2c"/2-2'; sinh277sinh(|) 

0:7(7^) = — 2e^"+^''sinh77sinh2?7 — e^''sinh477, 

y^{u) = 2e""^''sinh77sinh277-c''sinh477. (2.20) 

This 7?.-matrix is regular and unitary, with f{u) = xi{u)xi{—u). It is PT-symmetric and crossing- 
symmetric, with p = —677 — in and 

M = \ 1 I . (2.21) 



Diagonal solutions for K (u) have been obtained in [24]. It turns out that there arc three solutions 
without free parameters, being K^{u) = 1, K^{u) = F^ and K^{u) = F^ , with 






F^ = I g^{u) 1 , (2.22) 



where we have defined 



f (u) — cosh(-u — Sr/) ± i sinh(-?i), 9^{u) — cosh(-u + 3ry) =F i sinh(-u) (2.23) 

By the automorphism (2.7), three solutions K^{u) follow as K^{u) ~ M, K^{u) = 6*+ and K^{u) = G^ , 
with 

/ e"~4')/±(w) \ 

G± = /i±(u) , (2.24) 

V e-"+4V±(«) J 

where we have defined 

h^{u) — cosh(-u — Sr/) ^ isinh(-7i — Gr/)). (2.25) 

Here we will only consider three types of boundary solutions, one for each pair {K^{u), K^{u)) defined 
by the automorphism (2.7): (1,M), (F+,G+) and (F-,G-). 

The transfer matrix for the case (1,M) , whose corresponding spin-chain is f7q(su(2))-invariant [25] 
has been diagonalized by the analytical BA in Ref. [12]. For the other cases, (F^, G^) and (F^ , G^) ,the 
corresponding open chain Hamiltonians are not C/q(su(2))-invariant . Nevertheless, their transfer matrices 
were diagonalized by Yung and Batchelor in. [13] through the analytical BA with inhomogeneities. The 
energy eigenspectra for these three cases was also solved by Fireman et al [14] via the coordinate BA. 

Recently has been argued by Nepomechie [26] that the transfer matrices corresponding to these 
solutions also have the f7q(o(3)) symmetry, but with a nonstandard coproduct. 

It was in this model that Tarasov developed the algebraic BA for the three-state 19-vertcx models with 
periodic boundary conditions [17]. The analytical BA with periodic boundary conditions was presented 
in [6] and the corresponding coordinate BA was presented in [23]. 

3 Boundary Algebraic Bethe Ansatz 

In the previous section we have presented a common structure for the ZF and IK models. In this section 
we will turn to the eigenvalue problem for their double-row transfer matrix with integrable boundaries, 
named boundary algebraic Bethe Ansatz. 



3.1 The reference state 

The monodromy matrix T(u) (2.9) and its reflection T~^(—u) can be written as 3 by 3 matrices 

Tn{u) Tu{u) Ti3(u) \ / T.-^-u) T,'^\^u) T-,\~u) \ 

T{u) = I T2i{u) T22iu) T2z{u) , r-i(-M) = T^^\-u) T,',\-u) T,-^\~u) 

Tsiiu) T32iu) T,3{u) J \ T-,\-u) T~\-u) T.-^-u) J (3.1) 



where 



3 



T,a{u) = Y. 4^^("' ^) ® 4f.;'H": r;) ® . . . ® C^^lju, r;) (3.2) 



/ci,...,fciV-l = l 

(n) 

where £J ■ are 3x3 matrices acting on the nth site of the lattice, defined by 



(n) _ 







>C^i' =10 ye 



y5 
ye 





/OOyrX / ye \ 

4^) = , 4^) = 2/5 , 43^ = I X2 I . (3.3) 

\000/ \000/ 

Using the unitary relation in (2.2) we can sec that the reflected monodromy matrix T^^(—u) has the 
following matrix elements 

Ti/i-^) = TT^v E 411 (-"' -^) ^ 41(-«: -r?) ^ • • • C^Zlj-u, -ri). 

•^^"^ fei,...,fc„-i=i (3.4) 

Now we introduce the reference state 

|0)= ® ®---® I (3.5) 

The actions of T{u) and T^^(— u) on this state are 

/ x^{u)\0) * ** \ 

T{u)\0) = f^{u)T-\-u)\0)=\ x^{u)\0) *** (3.6) 

V x^iu)\0) J 

which give us, in the usual Bethe Ansatz language, the creation and annihilation operators for this 

reference state. Moreover, wc are dealing with reflection and in this case we have a double-row monodromy 

matrix deflned by 

Uiiiu) Ul2{u) Ul3{u) 

U{u)=T{u)K-{u)T-\-u)= { U2i{u) U22{u) U23{u) I (3.7) 

Usiiu) U32{u) U33{u) 



where K^ '{u) is a reflection matrix. 

For the diagonal case K^^'{u) ~ dia,g{kYi{u) , k22(u), k^^{u)), the matrix elements of U have the form 

3 

U,,{u) ^Y.^,,{u)k-,{u)T-,\-u), I, J = 1,2,3. (3.8) 

a=l 

It follows from (3.8) that we will need to know the commutation relations of the operators T{u) and 
T^^(— m) in order to get the action of U{u) on the reference state (3.5). Using (2.8) with u = —v we get 
the matrix relation 

T2\-u)Ri2{2u)Ti{u) = Ti{u)Ri2{2u)T2\-u) (3.9) 

Applying both sides of this relation on the reference state, we find the following relations between its 
matrix elements 



T,,{u)T^,\-u) |0) = h{u f' ^";^."^ ^"^ |0) (3.10) 



^13 y-<^) iw = j2\u)jjq-j |0) - /3(w)/i(«)-p^ |0) - fi{u)jW(^ 



ni{u)T,-,H-u) |0) = f2{u)^^^ |0) - f,{u)Mu)^^^ |0) - U{u)^^^ |0) (3.11) 



n,{u)T,-,\-u) |0) = Mu)^^^^j^^^^ |0) (3.12) 



where 



r, N 2/5(2u) ,, . Vi{2u) xi{2u)y^{2u) - x^{2u)yj{2u) 

Jim = — 77r\^ h{u) = — j—r, h[u) = — 7—- — 7—- , , . , 

xi[2u) xi[2u) xi[2u)Xi[2u) — X5[2u)y5[Zu) 

Mu) = x^{2u)yr{2u) - yl{2u) 

Xi(2u)x4(2u) — 2:5(2^)2/5 (2u) 

Using these relations we can get the action of all Uij{u) on the reference state 



Un{u)\Q) = fcnWl^lO) 






U22{u) |0) = h{u)U,,{u) |0) + [k^^iu) k^Mhin)] -jj^ 



C/33(«)|0) = [{h{n)- h{u)h{u))U^Au) + h{u)U22{u)m 

xfiu) 



+ [k-,{u) - k-,iu),Uu) - k-,{u)Mu)] ^^ |0) (3.14) 



and 



U,j{u)\0)=0, (*>j), U,j{u)\0)^{0,\0)}, {t<j) (3.15) 



In order to recover the usual BA structure we define news operators: 

Ci{u) = C/2i(m), V2iu) = U22{u) - h{u)Vi{u), B^{u) = C/23(m) 

C2(w) = t/3i(w), C3(u) = t/32(u), V^{u) = U:i3{u)-h{n}V^{u)-h(u)V2{u) 

and then the new double-row nionodroniy matrix has the form 



(3.16) 



Vi{u) Bi{u) B2iu) 
Uiu)^U{u)^\ Ci{u) V2{u) Bz{u) I (3.17) 

C2{u) Cz{u) V^{u) 



The action of U{u) on the reference state is now 



where 



Xi{u) |0) * ** 

W(u)|0)=( ^-2(^)10) *** I (3.18) 

^-3(^)10) 



x^^ {u) 
Xi{u) ^ Ki{u)-^ 



x^iu) 



X2{u) = [k^^{u)-k^^{u)fi{u)\ 



r(") 



xi^iu) 



Mu) = [k^,{u)-k^2{u)h{u)-k^Mhin)\^w^ (3.19) 

The transfer matrix tin) (2.10), with diagonal left reflection K''^^ = diag(fcj*j^, fc^.fc;^) has the form 
t{u) = k-^^iupiiiu) + k^2{u)U22{u) + k^^{u)Ui3{u) 

= ni{u)Vi{u) + n2{u)V2{u) + nz{u)v^{u) (3.20) 

where 

Q.i{u) = /c+(u) + /i(u)fc+(u) + /2(u)fc+3(u) 

^2(.u) = k+^iu) + f3{u)k^^{u) 

n^iu) = ^3+3(71) (3.21) 

It is now clear that for U{u) we have recovered the usual algebraic BA structure. Therefore we can 
look for states created by the action of the operators Bi (u) on the reference ^o which will be eigenstates 
of the transfer matrix (2.10). To do this we first recall the magnon number operator 

N 

M^J2^^k, Mfc = diag(0,l,2) (3.22) 



k=l 



10 



This is the analogue of the operator S^ used in the coordinate BA construction. The relation M'^m — 
m^fjn where m — N — S^, allows us to build states ^P^ such that t(u)^m = A„^„. Therefore, we can 
start the diagonalization of t(u) by considering all possible values of tti in a lattice with A^ sites. 
By the previous construction, ^Pq is our reference state |0 > , which is itself an eigenstate of t{u) 

i(u)*o = Ao(w)*o (3.23) 

with eigenvalue 

Ao{u) = [k+^(u) + Mu)k+,{u) + h{u)k+,(u)]K,(u)^j^ 



[fc+ (u) + /3(u)A:+ (u)] [fc22(u) - /cii(u)/i(u)] 






4" in) 
/^(«) 
This is the only state with m = 0. 



+k+,{u) [k^siu] - k^,{u)Mu) - k^,{u)Uiu)] %)^ (3.24) 



3.2 The one-particle state 

For TO = 1 we seek a state of the form 

*i(ui) = ^i(ui)|0). (3.25) 

Then the action of the transfer matrix t{u) on this state is 

t(w)*i(Mi) - niiu)Vi{u)Biiui) |0) + r!2(u)2?2(M)^i(Mi) |0) + n3iu)V3{u)Biiui) |0) . 

(3.26) 

Since we know the action of the operators 'Di{u) on the reference state |0), we need to arrange the operator 
products 

V^{u)Bi{ui), V2{u)Bi{ui) and V3{u)Bi{ui) (3.27) 

in a normal-ordered form [17] . 

We anticipate that, in general, the operator- valued function ^„('Ui,... ,u„) for a n-particle Bethe 
state will be composed by a set of normal-ordered monomials. A monomial is said to be in normal order 
if all elements Bi are on the left, and all elements Ci are on the right of the elements Vi. 

In order to get this normal ordering we recall that the double- row monodromy matrix U{u) satisfies 
the fundamental reflection equation 

TZi2{u - v)Ui{u)n2i{u + v)U2{v) - W2(w)7^l2(w + v)Ui{u)n2i{u - v), (3.28) 

where Ui{u) — U{u) ® 1, U2{u) — 1 ®U{u) and 7?,2i(m) = PTZi2{u)P. In the appendix we show how this 
equation (indeed a set of 81 equations for the three-state models) can be used to recast the non-normal 
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ordered operator products as the above into a linear combination of normal-ordered ones , which might 
be called fundamental set of generalized commutation relations, or shortly commutations relations. 

For the present case t{u)^! i{ui) can be computed with the aid of the following commutation relations 
(see the appendix) 

'Di{u)Bi{ui) = aii{u,ui)Bi{ui)'Di{u) + ai2{u,ui)Bi{u)Vi{ui) + ai3{u,ui)Bi{u)V2{ui) 

+ai4(M, ui)B2{u)Ci{ui) + ai5(u, ui)B2{u)C3{ui) + aie{u, ui)B2{ui)Ci{u) (3.29) 

T>2{u)Bi{ui) = a2i{u,ui)Bi{ui)T>2{u) + a22{u,ui)Bi{u)Viiui) + a23{u,ui)Biiu)D2{ui) 
+024 {u,ui)B3{u)T)i{ui) + a25{u,ui)B3{u)T>2iui) + a2e{u,ui)B2{u)Ci{ui) 
+a27iu,ui)B2iu)C3{ui) + a2siu, ui)B2iui)Ci{u) + a29iu,ui)B2iui)C3{u) (3.30) 

V3{u)Bi{ui) = a3i{u,ui)Bi{ui)T)3{u) + 032(1*, mi)Si(u)Pi(ui) + 033(1*, ui)Si(u)P2(ui) 
+034 {u,ui)B3{u)T>i{ui) + a35{u,ui)B3{u)T>2iui) + a3e{u,ui)B2{u)Ciiui) 
+a37{u,ui)B2{u)C3iui) + a3s{u,ui)B2{ui)Ci{u) + a3g{u,ui)B2iui)C3{u) (3.31) 

Introducing these relations in the eq (3.26) one gets 

i(u)*i(ui) = ni{u)Vi{u)Bi{ui) |0> + n2iu)V2iu)Biiui) |0> + n3iu)V3iu)Biiui) |0> 



= [aii{u,Ui)n.i{u)Xi{u) + a2i{u,ui)n2iu)X2{u) + a3i{u,ui)Q3{u)X3{u)] *i(ui) 
3 3 

+ [Xi{ui)^flj{u)aj2{u,ui) + X2{ui)^flj{u)aj3{u,ui)]Bi{u) |0) 

3 3 

+ [Xi{ui)^nj{u)aj4{u,ui) + X2{ui)^nj{u)aj5{u,ui)]B3{u)\0) (3.32) 

j=2 3=2 

So ^1^1(^1) will be an eigenstate of t(u) with eigenvalue 

3 

Ai(u,Mi) = y^Qj{u)Xj{u)aji{u,ui) (3.33) 

provided the following equations are satisfied 

XijUi) _ Ej=1^3Naj3(M,Ml) _ Yfj=2^3iu)(^J5{u,Ui) _ 
<^2(ui) Z]j=l%(")ai2(w,'"i) Z]j=2^j(")"i4(w, Ml) 

These are the Bethe equations for the one-particle state and the corresponding Bethe root ui does not 
depend on u. 
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3.3 The two-particle state 

For 171 — 2 there are two linearly independent states Bi{ui)Bi{u2) \0) and 62(^1) |0). Therefore we seek 
for eigenstates of t{u) in the form 

*2(mi, W2) = Bi{ui)Bi{u2) |0) + B2{ui)T{ui,u2) 10} (3.35) 

where r(wi,U2) is an operator- valued function. Next we will use the condition that ^2(^1, ^2) rnust be 
normal-ordered to find r('Ui,U2)- 

The first term in the right hand side of the eq (3.35) has its normal-ordcrd form given by the 
commutation relation: 

Bi{ui)Bi{u2) = Uj{ui,U2) [Bi{u2)Bi{ui) + GdJu2,Ui)B2{u2)Vi{ui) + Gd2{u2,Ui)B2{u2)V2{ui)] 

-GdAui,U2)B2{ui)Vi{u2) - Gd2{ui,U2)B2{ui)V2{u2) (3.36) 



where 



Xs,{Ui - U2)xi{ui ~ U2) - Xfi{ui - U2)V(>{ui - U2) 

U){ui,U2) = r ^ 

Xl{Ul - U2)X3(Ul - U2) 

Uj{u2,Ul)uj{ui,U2) = 1 (3.37) 

^ / N X6(m -U2) 3^2(2^2) .„ „„^ 

Gd^ [ui,U2) = T— T^r^ (3-38) 

X3[Ui - U2) Xi(2u2) 

GdA^uU2)^ — 7 ' r (3.39) 

X2(Ux +M2) 



Here we have used the following identities valid for both models, 

yej-u) ^ X3{u)xe{u) - X7{u)y(i{u) 
X3{-u) X3{u)x4{u) - xe{u)ye{u) 

and 

X2(2u) _ y5(u — v)x2{u + v) + X2{u — v)x<^{u + v)Ji{u) 
xi{2u) y^{u — v)xi{u + v) 

Now we can see that (3.35) is normal ordered if it satisfies the condition 



(3.40) 



(3.41) 



*2(m2, Ul) = ^(W2, 'Ui)*2(mi, W2) (3.42) 

This condition fixes r(Mi, U2) and, by construction, the unique candidate for the eigenstate of t{u) in the 
m = 2 case has the form 

*2(W1,W2) = B^{u^)Bi{u2) |0) + GdAui,U2)B2{u^)V^{u2) |0) + GdAui.U2)B2{u^)V2{u2) |0) 

(3.43) 
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The action of t{u) on this state reads 

t{u)-^2{ui,U2) ^ ni{u)Vi{u)Bi{ui)Bi{u2) |0) + n2{u)V2{u)Bi{ui)Bi{u2) |0) 
+n3{u)V3{u)Bi{ui)Bi{u2)\0)+Gd,{ui,U2)ni{u)Vi{u)B2{ui)Vi{u2)\0) 

+ Gd, {ui,U2)n2{u)V2{u)B2{ui)Vi{u2) |0) + Gd, {ui,U2)n3{u)V3{u)B2{ui)Vi{u2) |0) 
+ Gd,{ui,U2)ni{u)Vi{u)B2{ui)V2{u2)\0)+GdMuU2)^2{u)V2{u)B2{ui)Vi{u2)\0) 

+Gd,{ui,U2)n3{u)V3{u)B2{ui)V2{u2) |0) (3.44) 

In order to recast the right hand side of the above expression to its normal ordered form we will need the 
following commutation relations, in addition to those presented in the m = 1 case (see the appendix): 

Vi{u)B2{v) = bn{u, v)B2{v)Vi{u) + h2{u, v)B2{u)Vi{v) + h^iu, v)B2{u)V2{v) 

+bu{u,v)B2{u)V3{v) + bi5{u,v)Bi{u)Bi{v) + he{u,v)Bi{u)B3{v) (3.45) 

V2iu)B2iv) = 621 (m, v)B2iv)V2{u) + 622 (w, v)B2{u)Vi{v) + 623 (w, v)B2{u)V2{v) 
+624 {u,v)B2{u)V3{v) + 625 {u,v)Bi{u)Bi{v) + 626 {u,v)Bi {u)Bz {v) 
+b27{u,v)B3{u)Bi{v)+b2s{u,v)B3{u)B3{v) (3.46) 

V3{u)B2iv) = b3iiu,v)B2{v)V3{u) + b32{u,v)B2{u)Vi{v) + b33{u,v)B2{u)V2{v) 
+634 {u,v)B2{u)V3{v) + 635 {u,v)Bi{u)Bi{v) + 636 {u,v)Bi {u)B3 {v) 
+b37{u,v)B3{u)Bi{v)+b3s{u,v)B3{u)B3{v) (3.47) 

Ci{u)Bi{v) = cii{u,v)Bi{v)Ci{u) + ci2{u,v)Bi{v)C3{v) + ci3{u,v)Bi{u)C3iv) 

+ci4{u,v)B2{u)C3{v) + ci5{u,v)B2{v)C2iu) + cie{u,v)Vi{v)Vi{u) 
+cir{u, v)T>i{v)T>2{u) + cis{u, v)T)i{u)T>i{v) + cig{u, v)T)i{u)T)2iv) 
+cnoiu,v)V2{u)Vi{v) + cin{u,v)V2{u)V2iv) (3.48) 

C3{u)Bi{v) = C2iiu,v)Biiv)Ci{u) + C22{u,v)Bl{v)C3{v) + C23{u,v)Bi{u)C3{v) 

+C24 {u,v)B2{u)C3iv) + C25 {u,v)B2iv)C2 (u) + C26 {u, v)Vi{v)T)i{u) 
+C27{u, v)I)i{v)V2iu) + C2s{u, v)T>i{u)I)3{v) + C29(u, v)T>i{u)T>2{v) 

+ C21oiu,v)Vliu)V2{v) + C21l{u,v)V2{u)T>i{v) + C212{u,v)T)2{u)T)2iv) 

+C2i3iu,v)V3iu)Vi{v) + C2uiu,v)V3{u)Vi{v)) (3.49) 

After a straightforward calculation we obtain 

3 

t{u)'^2iui,U2) = [^flj{u)Xj{u)aji{u,Ui)aji{u,U2)]'^2iui,U2) 



i=i 



14 



+ [aii{ui,U2)Xi{ui)^rtj{u)aj2{u,ui) + a2i{ui,U2)X2{ui)^^j{u)aj3{u,ui)]Bi{u)Bi{u2) |0) 

3 3 

+ [aii{ui,U2)Xi{ui)^flj{u)aj4{u,ui) + a2i{ui,U2)X2{ui)^nj{u)aj5{u,ui)]B3{u)Bi{u2) |0) 

3=2 j=2 

3 3 

+ [an{u2,ui)Xi{u2)^flj{u)aj2{u,U2) + a2i{u2,ui)X2{u2)^^j{u)ajs{u,U2)]u!{ui,U2)Bi{u)Bi{ui) |0) 

3 3 

+ [aii{u2,Ui)Xi{u2)^iljiu)aj4{u,U2) + a2liu2,Ui)X2{u2)^^j{u)aj5{u,U2)]uj{ui,U2)B3{u)Bi{ui) |0) 
j=2 j=2 



+ [Xl{ui)Xi{u2)^flj{u)Hji{ui,U2) + Xi{ui)X2{u2)^^j{u)Hj3{ui,U2) 
3 3 

+X2{ui)Xi{u2)Y,n,{u)H,2{ui,U2) + X2{ui)X2{u2)Y,n,{u)H,4{ui,U2)]B2{u)\0) (3.50) 
where 

iJll(ui,U2) = ai4(u, Ml) {cie{ui,U2) +Cls{ui,U2)) +ai5(u, Ml) (c26(ui,M2) + C29(ui,U2)) 

+5i2(u, ui)Gdi (wi, U2) + tj(ui, u)aii(u, Mi)ai2(u, U2)Gdi {u, ui) 

Hi2{ui,U2) = aii{u,Ui) (C17(U1,U2) +Ciio(ui,U2)) +ai5(u,Mi) (C27(U1,U2) +C21l(ui,U2)) 

+5i3(m, ui)G<ii (wi, U2) + uj{ui,u)aii{u, ui)ai2{u, U2)Gd2 {u, ui) 

Hi3{ui,U2) = ai4(u, Mi)ci9(ui,U2) + ai5(w,'"l)c210(wi,'"2) + bi2{u,Ui)Gd2{ui,U2) 

+uj{ui,u)aii{u, ui)ai3 (u, U2)Gdi (u, ui) 

iJl4(ui,M2) = ai4(u, Ml)cill(ui,U2) + ai5(u, '"1)C212(U1,U2) +613(^,^1)6*^2(^1,1*2) 

+w(mi, u)aii(u, wi)ai3(u, U2)Gd2 {u, mi) (3.51) 
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and 

Hji{ui,U2) = aj6{u,ui) {cie{ui,U2) + cis{ui,U2)) + ajr{u,ui) (026(^1,^2) +229(^1,^2)) 

+bj2{u, ui)Gdi (mi, U2) + cj(mi, u)aji(u, Ui)aj2(u, U2)Gdi {u, ui) 

+aji{u, ui)aj4:{u, U2)di3{ui,u) 
Hj2{ui,U2) = aj(i{u,ui) {cn{ui,U2) + ciio(ui, U2)) + aj'j{u,ui) (027(^1, U2) + C2ii(mi,M2)) 

+6j3(m, ui)Gdi (mi, U2) + uj{ui,u)aji(u, ui)aj2{u, U2)Gd2 (u, ui) 

+aji{u, Ui)aj4{u, U2)di4{ui,u) 
Hj3{ui,U2) = aj6{u,ui)cig{ui,U2) + ajr{u,ui)c2ia{ui,U2) +bj2iu,ui)Gd2iui,U2) 

+uj{ui,u)aji{u, ui)aj3{u, U2)Gdi {u, ui) + aji (u, Ui)aj5{u, U2)di3{ui,u) 
Hji{ui,U2) = aj6(u,ui)ciii(wi,U2) + aj7(u, '"i)c2i2(ui, "2) + 6j3(u, ui)Gd2(ui, U2) 

+cj(ui,u)aji(u, Mi)aj3(u,M2)Gd2(u, Ml) + aji(u, ■ui)aj5(u, M2)di4(ui, u) (3.52) 

for i = 2, 3. 

Again, ^2(^1, U2) will be an eigenstate of t{u) with eigenvalue 

3 

A2(u,Mi,U2) = 2j^j(")'^j(")"jl('"''"l)"jl("'"2) (3.53) 

provided the following Bethe equations are satisfied 

-yi(m) ^/ a2i(Mi,M2) -^1(^2) ^. 021(^2,^1) .„ ^.^ 

^ e Ml ] T: , , ^9^2 ^ T 3.54) 

X2{Ul) aii(ui,U2) X2[U2) aii{U2,Ui) 

where <d{ui), i = 1,2 are given by (3.34). 

3.4 The three-particle state 

In this case we have to consider combinations and permutations of the states of the type BiBiBi |0) and 
B2B1 |0) . 

Let us consider the normal-ordered operator 

$3(1*1, M2,U3) = Bi{ui)<^2{u2,U3,) + B2{ui)^i{u2)Ti{ui,U2,Uz) + ^2 ("l )*1 (U3)r2 (wi , M2, "a) 

(3.55) 
with two exchange proprieties 

$3(W1, W2, m) = UJ{UI,U2)^3{U2, Ul, U3) == tj(u2, M3)$3(wi, "3, ""2), (3.56) 

which can be used, together the commutations relations, to find the operator valued functions 

ri(ui,U2,U3) = UJ{U2,U3) [aii(u3,U2)G<ii(ui,U3)^l(w3) + a2l(w3, ""2)0^2(1*1, U3)P2(W3)] 
r2(ui,U2,U3) = aii{u2,U3)Gd^{ui,U2)'Di{u2)+a2l{u2,U3)Gd2{ui,U2)'D2{u2) (3.57) 
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Note that 

Ti{ui,U2,U3) = UJ{U2,U3)T2{U1,U3,U2) (3.58) 

The eigenvalue problem for the 3-particle state is now reduced to the problem of finding the action 
of the operators !?„(«) , a = 1, 2, 3 on the state 

*3(mi, U2, U3) = <^3iui,U2,U3) |0) (3.59) 

To do this we will need recall the appendix to get more five commutation relations: 

'Di{u)B3{v) = Xii{u,v)B3{v)T>i{u) + xi2{u,v)Bi{v)T>i{u) + xi3(u,v)Bi{u)T>i{v) 
+xu{u,v)Bi{u)T>2{v) + xi5iu,v)Bi{u)V3{v) + xi(i{u,v)B2iu)Ci{v) 
+xir{u,v)B2{u)C3{v) + xis{u,v)B2iv)Ci{u) (3.60) 

'D2{u)B3{v) ^ X2i{u,v)B3{v)'D2{u) + X22{u,v)Bi{v)'D2{u) + X23{u,v)Bi{u)Vi{v) 
+2:24 (u, v)Biiu)V2iv) + X25{u, v)Bi{u)T>3{v) + X26{u, v)B3{u)T)i{v) 
+2:27 (u, w)^3 (u)2?2 (v) + a;28 {u,v)B3{u)T>3iv) + 2:29 {u,v)B2 {u)Ci (v) 
+X2iQ{u,v)B2{u)C3iv) + X2iiiu,v)B2iv)Ci{u) + X2i2iu,v)B2iv)C3iu) (3.61) 

'D3{u)B3{v) = X3i{u,v)B3{v)T>3{u) + X32{u,v)Bi{v)T>3{u) + X33{u,v)Bi{u)Vi{v) 
+2:34 (u, v)Bi{u)V2iv) + X35{u, v)Bi{u)T>3{v) + X3g{u, v)B3{u)T>i{v) 
+X37{u,v)B3iu)V2iv) + X3s{u,v)B3{u)T>3iv) + X3g{u,v)B2{u)Ci{v) 
+X3ia{u,v)B2{u)C3{v) + X3iiiu,v)B2iv)Ci{u) + X3i2{u,v)B2iv)C3{u) (3.62) 

Ci{u)B2{v) = yii{u,v)B2iv)Ci{u) + yi2{u,v)B2iv)C3{u) + yi3{u,v)B2{u)Ci{v) 

+yi4{u,v)B2{u)C3{v) + yi5{u,v)Bi{v)T)i{u) + yiG{u,v)Bi{v)T>2{u) 
+yir{u,v)B3{v)Vi{u) + yis{u,v)B3{v)V2iu) + yigiu,v)Biiu)Vi{v) 
+yiioiu,v)Bi{u)V2{v) + yiii{u,v)Bi{u)T)3{v) + yii2{u,v)B3{u)Vi{v) 
+yii3(u, v)B3iu)V2{v) + yiu{u, v)B3{u)V3iv) (3.63) 

^3(^)^2(1') = y2iiu,v)B2iv)Ci{u) + y22iu,v)B2iv)C3{u) + y23{u,v)B2{u)Ci{v) 

+y24:{u,v)B2{u)C3iv) + y25iu,v)Bi{v)I)i{u) + y26iu,v)Bi{v)T>2{u) 
+2/27 (w, v)Bi{v)V3{u) + y28{u, v)B3{v)Vi{u) + y29{u, v)B3iv)T>2{u) 
+y2ioiu,v)B3{v)T>3{u) + y2ii{u,v)Bi{u)T>i{v) + y2i2{u,v)Bi{u)V2{v) 
+y2i3iu,v)Biiu)V3{v) + y2i4:{u,v)B3{u)T)i{v) + y2i5iu,v)B3iu)V2{v) 
+y2i&{u,v)B3{u)V3{v) (3.64) 
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After a straightforward but (the reader should be advised) quite lengthy computation we obtain the 
following simplified expressions 

3 



-^Y\_Uj{Uj,Ui)[Xi{ui)aa2{u,Ui)Y\_ail{Ui,Uk) + Af2(Ui)aa3(w, Ui) J]^ a2l(ui,Ufe)];Bl(M)'I'2( 
i—1 j — 1 k^i k^i 



+ (1 - 5a.l) ^ W U}{Uj,Ui)[Xi{ui)aaA{u, Ui) Jl aii(ui, Uk) + X2{ui)aa^{u, Ui) JJ a2l{Ui, Uk)]B^{u)-^2{ui) 
i—1 j — 1 k^i k-^i 

2 3 

2—1 j'^-i+l 

-\-X2{u^)Xx{Uj)a2\(Ui,Um)a\\{Uj,U^)ila2{ui,Uj) + Xi{Ui)X2{Uj)aii{ui,Um)a2l{Uj,Um)Ha3{Ui,Uj) 

l~l j-1 

+ X2{Ut)X2{Uj)a2l{Ut,Um)a2liUj,Um)Ha4iUi,Uj)} JJ w(ufe,Ui) J]^ w(w; , Wj )S2 (u)*l (Um) 

fe=l ; = l5^i 

(3.65) 

for q; = 1, 2, 3 and jti ^ {*,.?}• In the above expression we have introduced the symbol Ui , which as usual 
means that the state ^2 has to be evaluated at spots different form Ui . 

By performing the action of transfer matrix t{u) on the state ^£'3(1*1, U2, U3), we get 

i(u)*3(ui,U2,W3) = i'^fla{u)Xa{u)Y]_aal{u,Ui)\ ■^3{ui,U2,U3) (3.66) 

\«=1 i=l / 

all the unwanted terms vanishing provided that the Bethe equations are satisfied: 

Xt(uj) „, , -i-r a2-\(Uj,Ui) , . , . 

_^^e(u.) n ■' (* = 1,2,3) (3.67) 

X2(Ui) .J-.-^ aii{u^,Uj) 

3.5 The n-particle state 

From the previous results one can seek for operator valued functions with a recurrence relation of the 
form 

$„(u,...,U„) = Si(ui)$„_i(u2,.. . ,Un) 

71 

+B2{ui)'^ Fi {Ui, . . . ,U„)$„_2(m2,--- ,Ut,... ,Un)'Di{ui) 



1=2 
i=2 



■B2(ui)^r2 (Ul,-- • ,U„)$„-2(U2,-- • ,U^,... , Un)T>2{Ui) (3.68) 



18 



It was shown in ([10]) that the above operator wiU be normal ordered satisfying n— 1 exchange conditions 

$„(W1, . . . , Ui, Ui+l, ... , U„) = Uj{Ui,Ui+l)^n{ui, ... , Ui+l, Ui, ... , u„) (3.69) 

provided that the functions Fa {ui, . . . , Un) are given by 

i—l n 

ri'^wi,.. . ,u„) = J|a)(Mj,Ui) Y]_ aai{ui,Uk)Gd^{ui,Ui), (a = 1,2) (3.70) 

j=2 k=2,k=ii 

Therefore the n-particle state will be given by 

*„(mi, . . . , u„) = *„(m, . . . , Un) |0) (3.71) 

and the action of the operators T>a{u) , a — 1, 2, 3 , on this state will be represented by 

n 
n i—1 n n 

-^^^Wuj{uj,Ui)[Xi{ui)ac2{u,Ui)^aii{ui,Uk) -^ X2{ui)ac,-:i{u,Ui)^a2i{ui,Uk)]Bi{u)'^ n-i{ui) 

i—l j—1 k^i k^i 

n i—l n n 

+ {l-5a.l)y^^\\ t^{U;j,Ui)[Xi(ui)aaA{u,Ui)\\ aii{ui,Uk)+X2{ui)aab{u,Ui)\\ a2l{Ui,Uk)]B'i{u^ 
i=l j = l k^i k^i 

n— 1 ^1 ^^ ^ 

+ ^ ^ lXi{ui)Xi{uj) J]^ aii{ui,Uk) W aii{uj,ui)Hai{ui,Uj) 

i=l j=i+l \^ k^i,j l=ii,j 

n n 

+ X2{ui)Xi{Uj) Y]^ a2l{Ui,Uk) Y\_ aii{Uj,Ui)Ha2{ui,Uj) 
k^i,j l^i,j 

n n 

+Xi{ui)X2{uj) Y\^ aii{ui,Uk) Y\_ a2i{uj,ui)Ha3{ui,Uj) 

k^i,j l^ij 

n n 

+ X2{ui)X2{Uj) Y]^ a2l{Ui,Uk) Y\_ a2l{Uj,Ui)Ha4:{ui,U 

x]^tj(Mfc,Ui) Y[ 0j{ui,Uj)B2{u)-^n-2{ui,Uj) (3.72) 

k=l l=l^i 

Finally, the corresponding n-particle eigenvalue problem will be 

/ 3 n \ 

t(u)*„(ui, . . . , U„) = ^ Q.a{u)Xa{u) J]^ a„l(M, Ui) *„(ui, . . . , u„) (3.73) 

\q=1 j=l / 

provided that the Bethc equations are satisfied: 

^ = e(..) n ^-^'^y ik-l,2,...,n) (3.74) 
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4 Explicit Solutions 

In this section explicit expressions of the eigenvalue problem are presented for both models. First we recall 
the appendix to get the coefficients aij{u,v) which appear effectively in the Bethe Ansatz expressions 
(3.73) and (3.74): 

Xl{v — u) X2{u + v) 



aii{u,v) 



a2i{u,v) = io{u,v)[ 



X2{v — u) Xl(u + v) 

xi{u + v)x4(u + v) — 2:5 (m + v)y^{u + v) 



Xi{u + v)x2{u + v) 



, ^ X2{u - v) X2{u + vf - Xfi{u + v)yQ{u + v) 
asi u,w) = — — ■ — ^ — - — ■ — 4.1 

X^{U-V) X2[U + V)X^{U + V) 

For the factor with the boundary contributions Q{ui), we will consider only the simplest expression 

Q. . _ 1^2 (u)a25 («,•»») + ^3 (M)a35(M,'»i) ,.^. 

Q.2{u)a2i{u,Ui) + Vlz{u)a^i{u,Ui) 
where the coefficients are given by 

xq{u ~ v) xz{u + v) xi:i{u + v) 

"^^("'") - W-.).2(u+.) '^^^"^^l(^rT^J 

a..M ^ -^i4^ (4.3) 

X2\U + V) 



034 (w,w) = !z{u)[fi{v 
+h{v) 



X(i{u + v) xq{u - v) X3{u + v) 
X2{u + v) X3{u — v) X2{u + v) 

ye(u - v) xq{u + v)ye{u + v) - x^iu + v) 



x^{u — v) 2:2 (u + w)x3 {u + v) 

xeju - v)ye{u - v) - xlju - v) yeju + v) 

xl{u — v) X2{u + v) 

xe{u + v) yaiu - v) xe{u + v)ye{u + v) - xliu + v) 
a35[u,v) = /3 M — -, — ■ — r^ -. -. — ■ — ^ — -, — ■ — ^ 4.4 

X2(U + V) X3[U-V) X2[U + V)X3{U + V) 

4.1 ZF Model 

Substituting the matrix elements of the TZ matrix and of the K matrices for this model we get the 
following expressions for the terms with boundary contributions 



X2iu) - 



(3 sinh M + 2 cosh u x{ (u) 
/3 sinh u — 2 cosh u p^{u) 
sinh 2u [i sinh(u + 2?7) - 2 cosh('u + 2-q) x\^ (u) 
sinh(2u + 277) /3 sinh u — 2 cosh u P^(u) 

sinh(2u — rj) P sinh('u + 77) — 2 cosh(w -\- rf) (3 sinh(u + 2rj) — 2 cosh(M + 2rf) x\^ {u) 
sinh(2M + r]) P sinh {u — rj) + 2 cosh(w — rj) (3 sinh u — 2 cosh u p^{u) ^ 

(4.5) 
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Oi(u) 



sinh(2w, + 3?7) a sinh m — 2 cosh u a sinh(M + 77) — 2 cosh(M + 77) 

sinh(2u + rj) a sinh(u + 2ri) — 2 cosh(u + 2ri) a sinh(M + 77) + 2 cosh(M + 77) 
sinh(27i + 277) a sinh 7i — 2 cosh M 

sinh 2u a sinh(u + 277) — 2 cosh(u + 277) ' 



and 



, , a sinh M + 2 cosh u tA r\ 

^^"^ " " a sinh(u + 277) - 2 cosh(u + 277) ' ^ ' 



sinh(2ui + 27y) a sinh(7ii + rj) + 2 cos]i{ui + 77) 
sinh(2'Ui) a sinh('Ui + 7;) + 2 cosh(ui + 77) 



The coefhcients aii{u, v) are given by 

sinh(u + v) sinh(u — v — 2r]) 



aii{u,v) 
a2i{u,v) 
031 (u,w) 



sinh(M + V + 277) sinh(7i — v) 

sinh(u + v) sinh(u + v + 877) sinh('u — 77 + 77) sinh(u — v — 2rf) 
sinh(u + 7; + 277) sinh(u + f + 77) sinh(u — 77 — 77) sinh(u — 7;) 
sinh(7/ -77 + 77) sinh(M + 7; + 877) 



sinh(7/ -77 — 77) sinh(7i + 77 + 7;) 
021 (u, 7;) sinh(M + 77 + 87;) sinh(u -77 + 77) 



(4.8) 



ail (77, v) sinh(77 + f + 77) sinh(77 — 77 — 77) 

It means that the 77-particle state ^„({7ii}) (8.71) is an cigcnfuction of the ZF transfer matrix with the 
eigenvalue 

sinh(7i + Ui) sinh(7i — u^ — 277) 



A„(7i,{77j) = Xi{u)^i{u)W- 



-A'2(77)r!2(7i)J] 



^_^ sinh(7i + 77i + 2?;) sinh(7i — Ui) 

sinh(77 + Ui) sinh(77 + 7^^ + 877) 



X 



n 



'-J- sinh(77 + 7ii + 277) sinh(77 + 77^ + 77) 
sinh(77 -7^2+77) sinh(7i — Ui — 2r]) 



*- sinh(77 — Ui — rj) sinh(77 — 77^) 

-^ Sinh(77 -77^+7?) Sinh(77 + 7^^ + 87?) 

+^"1 7i fii 77 1 1 ^— T^TT — \ — ; — r ^^■^> 

^j'- sinh(77 — Ui — rj) sinh(7i + Ui + rj) 
provided that the parameters Ui satisfy the Bethe equations 

sinh(7ii + 277) \ a sinh(7ii + 77) + 2 cosh(?7i + ??)/? sinh(7ii + 277) — 2 cosh(Mi + 2ri) 



sinh(7ii) J a sinh(7ii + 77) — 2 cosh(77i + r;) /3 sinh 7ii + 2 cosh 77^ 

-p-r sinh(7ii + TXj + 877) sinh(7ii — Uj + 7;) 
-'■-'- , sinh(77i + Uj + rj) sinh(77.i — Uj — rj) 

i = 1,2, ...,77 (4.10) 

Putting these expressions in a symmetric for (7^^ -^ Ui—rj) one can sec that we have recovered the previous 
resuhs obtained by the fusion procedure [22, 13]. 
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4.2 IK Model 

In this model the K matrices have no free parameters and we have to consider two cases: the quantum 
group invariant and the non-quantum group invariant. For both cases the aii{u, v) coefficients are given 

by 

aii{u,v) 



sinh(i(u + i;)) sinh(^(M — w) + 2?7) 



smh{^{u + v) — 2rj) sinh(i(u — f)) 



a2i(u,w) 



sinh(|(u + v) — 477) sinh(i(?i — v) — 2vi) cosh(i(M — v) + -q) cosh(i(u + w) — 77) 
sinh(i(u + f ) — 2ri) sinh(i('u — v)) cosh(i(u — v) — rf) cosh(i(u + f ) — irj) ' 

031 (m,w) 



cosh(i(M — w) — Sr/) cosh(i(u + v) — 677) 
cosh(i('u — v) — f]) cosh(i(u + v) — 877) ' 



a2i{u,v) sinh(i(u + v) — 4t]) sinh(i(u — v) — 277) cosh(i(M — u) + 77) cosh(i(u + v) — rj) 
aii{u,v) sinh(i(M + w)) sinh(i(u - i;) + 277) cosh(i(M - u) - 77) cosh(i(u + u) — St;) (ah) 



Though, the boundary contributions arc different for each case: 
4.2.1 The quantum group invariant case 

In this case we have 

X.(u) = 2, sinhz. x^ju) 
^^ ' sinh(u-2r/) p^ {u) 



X,{u) = e'\J^^'' .J^^^rlt"\<"} (4-12) 



sinh('u — Arf) cosh('u — irf) p (u) 

sinh('u — 677) cosh('u — 77) 
sinh(w — 277) cosh(7i — 877) 
-2,; sinh('u - 677) 



sinh(u — 477) 
nsiu) = e-^" (4.13) 



and 



eM-c-^'' ""^^":-^^^ ..1,2,...,7^ (4.14) 

smh Ui 

Therefore, the 7i-particle state ^n({Mi}) is an eigenfuction of the IK transfer matrix t(u) with eigenvalue 
, r i\ sinh(u — 67;) cosh(u — 7;) x?^(u) -p-r , , 



sinh(u — 2?7) cosh(M — 877) p^{u) 

sinh(u — 677) sinhu x'^'^ (u) -n , , 
sinh(7i — 477) smh(M — 2?7) p" [u) ^j- 

sinhu cosh('u — 577) cco^ (m) -r-r , , ,,^^s 

smh(u — 477) cosh(M — 877) p" (u) -»--•- 



i=\ 
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provided that the parameters Ui are solutions of the Bethe equations 

sinh(^Ui — 277) \ -A- sinh(i(?ii + Uj) — Aif) sinh(i(uj — Uj) — 2r]) 

sinh(i?ii) / Ji j,i sinh(i(Mj + Uj)) sinh(i(?ij - Mj) + 277) 

cosh(i(ui — Uj) + ry) cosh(i(ui + Uj) — 77) 
cosh(i(ui — Uj) — rj) cosh(i(ui + Uj) — 877) 
i = l,2,...,n (4.16) 

Putting these relation in a symmetric form (ui -^ Ui + 277) one can see that we have the results obtained 
by ([12, 13]) using the analytical Bethe Ansatz. 

4.2.2 The non-quantum group invariant cases 

These were the cases considered by Guang-Liang Li, et al [10], where a spurious dependence on the 
spectral parameter could exist. Here we shall see that this dependence of the Bethe equations on the 
spectral parameter is inexistent.. 

As the correspondent ii'-matrix solutions are complex conjugated, we will consider here only the case 
of iF+,G+), defined in Section 2. 

A way to turn more direct the comparison of our results with those previously known is to borrow 
the notation presented in [13]: first, the contributions from the boundaries have the form 



X,{u)n,{u) = a{u)^j^, X,{u)n,{u)^f3{u)^j^, X,iu)n^{u) = j{u)^L^ 
p" [u) p" (u) p" (u) 



(4.17) 



where 



and 



a(") = ^'^ y! e+(M)sinh(M-677) 

Jj\ \^ Li j 

sinhti cosh?/ — 7 sinh(7i — 277) 

piu) — a{u) . -—- r-r-, ;:rT 

sinh(7/ — 477) cosh 77 + t smh(?i — 277) 

, , , s£~(u) sinhTi smh(u ~ 2ri) 

7(7i) = a{u)^^-— -^-^, ^ (4.18) 

4^(77) sinh(7i — 477) sinh(7/ — 077) 



C=^(7i) = 2 cosh(7i - 377) ±2ismh2r) (4.19) 



Second, the boundary factors 0(77^) are given by 

_ e~"' sinh(7/t - 2ri) cosh(|(77 + Uj) - 377) + ismh{^{u - Uj) - 2r]) (Aon\ 

sinhui cosh(i(7i — Tii) — 77) + 7sinh(i(7/ + Tii) — 477) 

Here we note that (4.20) is the expression obtained in [10] for their /3{u, Ui) factors. From this result it 
follows the dependence on u in their Bethe equations. However, analyzing the (4.20) expression carefully, 
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more precisely, using the identity 

'^n) -A- i sinhf if?;, — 7/,:^ — 9n\ rnsh n — i «inln ('■?/. — 9ti^ 

(4.21) 



cosh(i(u + Ui) - 3ri) + i sinh(i(u - Ui) ~ 2r]) cosh?] - i sinh(uj — 2ri) 



cosh(i(u — Ui) — f]) + isinh(i(M + Ui) — 4rj) cosh(Mi — Sry) 

one we can see that this u dependency is only apparent. Therefore the correct expressions for 0{ui) are 
given by 

e(..) ^ c-"- ^^"^("J - y '^"^^^ -^su,Hu. - 2rj) ^^^^^^^^ ^^^^^ 

sinh(Mj) cosh(ui — Srj) 

Next, to recover our notation we will need of the following identity 

^2(m) ^q/x smhu cosh(M-377) 
fli(u) sinh(u — 477) cosh(u — 77) 

Therefore the n-particle state '^n{{ui}) is an eigenfunction of the IK transfer matrix with the boundary 
contribution (F+ , G+ ) , and has eigenvalue 

A„(u, {ui\) = a{u) '^ TT aii{u, m) + I3{u) ^ TT a2i{u, m) + 7(^) V^ \ IT «3i(u, Ui) 

P^N f^^ P^(«) t\ Z' (") t=i (4.24) 

provided that the Ui parameters satisfy the Bethe equations 

sinh(iuj — 2ry) \ (i{ui) sinh(7ii - 477) cosh(7ii - 77) 



sinh(iui) / a{ui) sinhtii cosh(ui — Sr;) 

-A- sinh(i(7ii + Uj) - Ari) sinh(i(7ii — Uj) — 2t]) 
j^ij^z sinhi^im + Uj)) sinh(i(M, -7ij) + 27/) 
cosh(i(7ii — Uj) + rj) cosh(i(ui + Uj) — rj) 



(4.25) 



cosh(2(ui — Uj) — rj) cosh(2(ui + Uj) — Srj) 

Writing the above Bethe equations in a symmetric form (ui -^ Ui — rj) one can see that is exactly the 
result obtained in [13] using the analytical BA. 

5 Conclusion 

The main result of this paper may be summarized by saying that the Sklyanin algebraic Bethe Ansatz 
is unique and consistent for the ZF and the IK vertex models with boundaries, corroborating all previous 
investigations in these models. Thus the Bethe vectors here obtained can be used for instance to obtain 
the eingenvalue problem for the corresponding Gaudin models [27] with boundaries and their associated 
systems of Knizhnik-Zamolodchikov equations [28]. Also, our results encourage the believe that other 
19-vertex models , such as the osp(2|l) model and the s^(2|l)(^' model in the presence of boundaries can 
also be solved by the algebraic Bethe Ansatz. In fact, preliminary results [29] indicate that this is indeed 
true. 
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It was in the IK model that Tarasov developed the algebraic BA for the three-state models with periodic 
boundary conditions [17]. Here, with the aid of previous works [16, 10], two of the three-state 19-vertex 
models have their algebraic Bethe Ansatz derived using a generalization of the Tarasov's approach. 
The algebraic BA for n-state models with periodic boundary conditions was developed by Martins in 
[30]. Therefore we believe that the Martins's approach can be generalized to include the diagonal open 
boundary conditions. 

Acknowledgment: This work was supported in part by Fundagao de Amparo a Pesquisa do Estado 
de Sao Paulo-FAPESP-Brasil and by Conselho Nacional de Desenvolvimento-CNPq-Brasil. 
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A The Commutation Relations 

The equation (3.26) gives us the commutation relations for the matrix elements of the double-row mon- 
odromy matrix which play a fundamental role in the algebraic Bcthe Ansatz. Here we present the 
commutations relations and their coefficientes using a compact notation. Recall that the operator 

/ Pi(«) Bi{u) B2(u) \ 
U{u) = Ci{u) V2{u) Bsiu) , (A.l) 

V C2{u) C3{u) Vsiu) J 

satisfy the fundamental reflection equation 

ni2iu - v)Uiiu)n2i{u + v)U2{v) =U2{v)'Ri2{u + v)Ui{u)'R2i{u-v). (A.2) 

Substituting (2.14) and (A.l) into (A.2) we get 81 equations involving products of two matrix elements 
oihi{u). These equations can be manipulated in ordem to put the product of pairs of operators in the 
normal ordered form. To do this we shall proceed in the following way. First we denote by E[i,i] = the 
(i,j) component of the matrix equation (A.2) and collect them in blocks B[i,i],i — l,...,'b,i — i,...,10 — i, 
defined by 

B[i,j] = {F,, = E[i,j], /,, - E[j, i], FF,, - E[10 - i, 10 - j], //,, = E[IQ - j, 10 - i]} 

(A.3) 

>From these blocks we can see that the pair {Fij , /^ ) as well as [FFij , f fij) can be solved simultaneously. 
We introduce the notation 

D,^V,{u), d, = V,{v), B,^B,{u), h^B,{v), C,^C,{u), c,^C,{v) (A.4) 

for the operators of the double-row monodromy matrix and 

Xi^ Xi{u + v), Yi^yi{u + v), Xi ^ Xi{u-v), y-i = Viiu-v), (A. 5) 

for the Boltzmann weigths and 

{Z}ij = {z}ij{v,u), {z}y = {z}y(u,w), F^ = f^{u), fi = fi{v). (A. 6) 

for the coefficientes of the commutations relations, where 

{Z} = A,B,C,D,E,X,Y and {z} ^ a,b,c,d,e,x,y. (A.7) 

Taking into account these simplifications, we will indicate the pair {Fij, fij) or {FFij, f fij) for which 
the corresponding normal ordered relations were obtained: 
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• (i^l4,/l4) 

Dibi = aiibiDi + ai2Bidi + ai3Bid2 + ai4B2Ci + ai5B2C3 + aiQb2Ci 
Cidi = AndiC'i + A12D1C1 + A13D2C1 + AiiBiC2 + A15B3C2 + AiebiC2 (A.8) 

where the coefficients are 

a;2 -^1 -^1 2:2 Xi Xi 

. _ VbXe _ Xr xiXe 

Au - - — T^, Ai5 - --— , A16 - — -^. (A. 9) 

a;2 Xi Xi X2 Xi 



(FnJ: 



17; 



Dib2 = 61162-Di +5i2-B2c'i +^13-82^2 + ^14-82^3 + ^15-8161 +616-B1&3 

C2di = BndiC2 + B12D1C2 + B13D2C2 + B14D3C2 + Bi^C'iCi + B1QC3C1 (A.IO) 



where 



J, ^1 ^3 J, 17 2/6 ^6 j^ ^7 2/7 ^3 p ye -''^e p -'^7 

-Dli = — -— , -012 - -i^l — -r; r2— -r^, -B13 — -^3-r^, 

X3 Xi X3 Xi Xi X3 Xi X3 Xi Xi 

Sl4 = -f, ^- = -^f, ^---f- (A-11) 

Xi 2:3 Ai Ai 

• {FF3eJf36) 

Dib3 = xiib3Di + xi2biDi + xi3Bidi + xi4Bid2 + xi5Bid3 + X1QB2C1 

+X17S2C3 + a;i8fo2C'i 
Csdi = XiidiCg + Xi2rfiCi + X13D1C1 + X14D2C1 + X15D3C1 + XieBiC2 

+Xi7B3C2+XisbiC2 (A.12) 

with the foUowing coefficients 

_ X2 X3 2/5 Yq ye Y5 

-^11 — — -r^, -^12 - — -r; T^^ii; 

2:3 A2 2:3 A2 2:3 A2 

V _ p 2/6 ^4 rp Xe yr Yq ye Y5 

X3 X2 X2 X3 X2 X3 X2 

V j^ Xq ye X4 ye Y5 
-^14 — -^3-r^ -r; T^^i3, 

A2 2^3 X2 X3 X2 

V Xe y7 ye Y5 

X2 X3 X3 X2 

X,7 ^ -^f -^fA,„ X,,.^-^fA,e. (A.13) 

X3 X2 X3 X2 X3 X3 X2 

Note that for each pair of equations the corresponding commutation relations are related by inter- 
changing 

M <^ w, Di ^ di, Bi'r^a, Ci ^ bi (A. 14) 
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{F24, /24) 

Cibi = ciibiCi + ci2biC3 + ci^Bic^ + CUB3C3 + ci5b2C2 + ciediDi + cndiD2 

+cisDidi + cigDid2 + ciioD2di + ciiiD2d2 (A. 15) 

with the foUowing coefRcients 

X4 2:5-^6 J/5 ^6 ^7 ^5 

Cll = -T^, Ci2 = -r^, Ci3 = -— , Ci4 = --r^, Ci5 = -— , 

Ai X2 Xi X2 Xi Xi Xi 

Ai 2:2 A 1 2:2 Ai 2:2 A 1 Xi 

Cl9 = -( — 1^ + ^11^)' Clio == -/l^r^, Ciii = --— . (A.16) 

X2 Xi Xi Xi Xi 

The pairs (Fi6,/i6), {FisJis), (i^i^ie, //le), {FFi8,.ffi8) and F19 form a closed set of equations 

from which we have derived the following commutation relations 

-B261 = 6116152 + 61262-81+61362-83, 
C1C2 = -B11C2C1 + -E12C1C2 + -B13C3C2 

-B3&2 = 64162-83 + 6426ii?2 + 6436352 

C2C3 = -B41C3C2 + -E42C2C1 + -B43C2C3 

-8162 = 62162-81 + 62262-83 + 6236i_B2 + 62463-82 

C2C1 — E21CIC2 + E22C3C2 + E23C2CI + E24C2C3 
-8263 = 63163-82 + 63261B2 + 63362B1 + 63462-83 

C3C2 = -E^3lC2C'3 + -E'32C2C'i + -B33C1C2 + -B34C3C2 

S262 = 62B2, C2C2=C2C2 (A.17) 



where 



2:2^2 ys 2:2X6 

611 = — T^, 612 = — , 613 = —^r^- (A.18) 

2:1 A3 Xi Xi A3 

j;, 2:2 X2 2^2 >6 j^ X5 / A 1 n\ 

2^1 A3 Xi A3 Xi 

_ X2 X2X3 X5{xj + xl ~ xr^y^) X2XQ 

621 — -TTt T7~7^i 622 



2:1 X| - XgFe ' a;i(x| - 2:52/5) X| - XeYf 



6 



2:1X5 ^6^6 , 2:5 Xj 0:12:2 X3X6 
"^23 — -^ T72 y ,^ H -^ , 624 — 2 "^2 v v l^.^Uj 

2^2 - 2:5^5 A| - ATeFe 2:1 X^ - XqYq x^ - x^y^ A| - XqYq 

Xi2:2 X2X3 _ x^{x\ + xl ~ x^y^) X2YQ 

C31 — -^ -rp) TrTTi C32 



xl - X5J/5 X| - Xei^e ' ' 2:1 (x^ - X5y5) X| - XaYa 

X2 X^Yq X5 XgYe 2:1X5 X| /A r,1^ 

Xl X| - ATeFe 2:1 X^ - X^Ya x^ - x^y^ X^ - XeYe 
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(-F45,/45) 



I?2&i = a2ibiD2 + a22B\di + a23Bid2 + a^iB^di + a2bBs,d2 + a2&B2Ci 

+027-6223 + a28^2C'l + 029^2^3 
Cid2 == A2ld2Ci+A22DiCi+A2^D2Ci+A2iDiC^+A25D2Ci+A2QBiC2 

+A27B3C2 + A2sbiC2 + A2963C2 (A.22) 



with 



_ X4 X4 xe Xe X4 Y5 , ^e ^ 

^21 — — -r^ + -fs T^ -\ T^^13 H ^14 

X2 X2 X2 X2 X2 X2 X2 

0-22 = -Ji — -T? {Fi-\ -r^)ai2 H -r^An-l A12 

X2 X2 X2 X2 X2 X2 X2 

ys ^4 . i^ ^5 ^5 ^ ^ ^6 , a;6 Xe 

a23 = t; (-^1 H T^)ai3, a24 = -,/it^ H -''-ii, 025 = -"t^ 

X2 X2 X2 X2 X2 X2 X2 

a26 = -(-Ti H -— jai4 H -— Aig H Ais 

2^2 X2 X2 X2 X2 

Vb X5 / p , y5 >5 ^ ye . J. , ys Y5 X4Y5 ^^ 

a27 = t; [Fi^ -—)air,, 028 = (-Fi H -r^jaie H -t^Am H Aie 

X2 X2 X2 X2 X2 X2 X2 X2 Xq X2 

X4X5 X4Y5 . ^xe 
0.29 = — -rr -^ -rr^i5-\ A 17 (A. 23) 

X2 Xq X2 Xq X2 



(-P'lSi /15) 



Bibi = eoibiBi + 602^2-02 + 603^2-01 + 604-82^1 + eo5B2d2 

C\ci = EoiCiCi + Eo2d2C2 + Eo3diC2 + E04D1C2 + E05D2C2 (A.24) 



with 



X3X4 - xeve _ X3X4 - xeye Xq 
eoi — , eo2 — T^ 

XiXs X1X3 X2 

xsVe - xeyr X3 X3X4 - x^ye Xe 

X1X3 X2 XiXs A2 

X6 X3 Xq Xq 

eo4 = — -^-/il^' eo5 = --r^ (A.25) 

2:3 A2 A2 A2 

These conimutaion relation play a special role in the construction of the n-particle states, as men- 
tioned above. 

{FF2e,ff 26) 

S1&3 = diib:iBi + di2biBi + di3b2Di + ^1462-02 + ^15-82^1 + dieB2d2 + dnB2d3 
C3C1 = D11C1C3 + D12C1C1 + Di3diC2 + Di4d2C2 + D15D1C2 + D1QD2C2 + DnD3C2 

(A.26) 
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where 



and 



where 



n Xg {ysYe + X2Y5Bi5)Eni - ysYe 



i?13 
i?15 



Xi + Y^B,e' X2(X4 + nSi6) 

(j/s^e + 3:2^^15)^03 - fiy5X2 + X2Y5Bn 

X2{X^ + Y5Bie) 

(2/5^6 + X2Yr,Bi5)En2 - ybX2 

X2{Xi + Y5Bifi) 
[ybYa + 3:2^5^15)^04 + Fiy^X2 + ^2^:2^5 + X2Y^Bi2 
X2{X^ + Y5Bie) 

{ysYe + X2Y5Bi5)Eo5 + t/5^2 + ^^33^2^5 + 3:2^5^13 

X2(X4 + Y5Bie) 



"" ^ -t^ <--> 



-63^1 = ^21^1-83 +d22&l-Bl +d23&2-Dl +^24^2-02 +^25^2-03 + (i26-S2(il + (^27-82^2 
C1C3 = D21C3C1 +-D22C1C1 +D23C?1 (^2 +-024^2^2 + -D25cf3C'2 + -D26-D1C2 + -D27-D2C2 

(A.28) 

1^21 
C^23 
fi24 
d26 

d27 = -^^-^^eo5 (A.29) 



X4 Y5 
X3 X3 


■Bie, a 


hi 


y^Ye 
X2X3 


y^Ye 

V ^01 + 
3:2^3 


Y -^15 
-^3 




^^2X3 




- 


y^ Ye 

V ^03 
3:2 ^3 


+ ^^ B,2 






y5 ^2 ^ p ^5 
3:2 X3 A3 


ys 
3:2 


^6 ^ 

•^^eo2 + 


^ Bi3, ^25 
X3 


^t- 


X3^- 


. 2/5 X2 

/l V 

3:2^3 


y5 Ye 

X2 X3 


eo4 


X3 






y5 X2 

X2X3 


y,Ye 
X2X3 
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{F28, fis) 



_D2&2 = &2I&2-D2 + 622-82^1 + 623-62^2 + b2iB2d3 + b25Bibi + 626^1 &3 

+627S3&I + b28B3b3 
(72^2 = B2ld2C2 + B22D1C2 + B23D2C2 + B24D3C2 + B25C1C1 + B26C3C1 

+B27C1C3 + B28C3C3 (A.30) 
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where 



„ ^ , X5 X3 Ye J- D I ^5 "^3 J I ^6 

-D21 = -L H -T^a27 + -^^05, -D22 — — /1-D12 H -T^a23 + ——eo3 

X2 X2 X2 X2 X2 X2 

B23 = —flBi3-\ -T^d24 + -^eo2, -B24 = — /1-B14 H '--r^d25 

X2 X2 X2 X2 X2 

B25 = -.fiBi5-\ -r^d22 + -r^eoi, -B26 = -/1-B16 H -r^<^2i 

X2 X2 X2 X2 X2 

B2, = -^^, B2, = -^ (A.31) 

X2 X2 X2 

C361 = C2161C1 +C22&1C3 +C23-B1C3 +C24S3C3 +C25fo2C'2 +C26dl-Dl 

+C27diD2 + C2sdlD3 + C2gDidi + C21oDid2 + C2llD2di + 0212-02^2 
+C213-D3rfl + C214D3d2 

Cib^ — C2161C1 + C22fe3C'i + C23B3C1 + C24B3C3 + (72562(^2 + C2ediDi 
+C27d2Di + C2sd3Di + C29Didi + C2wD2di + C2iiDid2 

+ C213Did3 + C2uD2d3 (A.32) 

where 



C21 = 


ye Y^ X4 ye Y^ 

= T^(cii - 1), C22 = -r^ci2, C23 = 

X3 X2 X3 X3 X2 


_yT_ _ 
X3 


ye Y^ 
2:3 X2 


C24 = 


ye X5 ye Y5 ye Xi ye Y5 

V V "^^4, C25 Ci5 
X3 X2 X3 X2 X3 X2 X3 X2 






C26 = 


2:4 le , rp xeXz ye Yj ye Y5 
X3 X2 X3 X2 X3 X2 X3 X2 






C27 = 


^ xe X3 X4 Ye ye I5 xe X3 

- ^3 ^ + ^ V Ci7, C28 = „ 

a;3 X2 0:3 ^2 X3 X2 X3 X2 






C29 = 


f yi Ye yeXi f „ Xe ye >5 
h y hi'i- ^ ./i^2 „ „ C18 

•^3 -^2 2:3 X2 X2 X3 X2 






C210 = 


y? Ye ^ yeXi Xe ye Y^ 
X3 X2 X3 X2 X2 X3 X2 






C2II = 


r y&Xi r j^ Xe ye Y5 
X3 X2 X2 X3 X2 






C212 = 


ye X4 Xe ye Y;^ Xq 

Y J"i Y Y "^iii' ^^213 - ./i V ' 

X3 X2 X2 X3 X2 X2 


C214 = 


Xe 

X2 


• {F F23, If 23) 









-^ (A.33) 



-D361 = a3ibiD3 + a32Bidi + a33Bid2 + az^Bzdi + a3^B3d2 + a3eB2Ci 

+037-6223 + a2s.h2Ci + a29b2C3 
Cid3 = Asidsd + A32D1C1 + A33D2C1 + A3iDiC3 + A35D2C3 + A3eBiC2 

+A37B3C2 + A2sbiC2 + A29b3C2 (A.34) 
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where 



_ X2 X2 X2 Ye 
X3 X3 X3 A3 

f yr Y5 ys ¥7 X2Yfi 
az2 = -Ji — T^ H T^^ii - ((^1)022 - ((^2)012 H -r^Ai2 

2:3 A3 2:3 A3 X'i A3 

033 = -^17 (Q2)ai3 - (Qi)a28 

2^3 -^3 

„ 2/6 ^2 ,^ . , 2^2 >6 ^ ye X2 , , 

034 = -/I — t; (Qi)a24H T^-^ii, 035 = -^ (^1)025 

2:3 A3 xa A3 2:3 A3 

y^ ^7 - ,^ . f^ . X2Ye 

2:3 A3 2:3 A3 

037 = -^-r? (Q2)ai5 - (Qi)a27 

2^3 A3 

«38 = — -rr^u - [Q2)ai6 - [Ql)a28 ^ T^^16 

2^3 A3 2:3 A3 

a39 = — T^H T^Ai5 - (Qi)a29 H -rrXn (A. 35) 

Xz A3 Xs A3 2:3 A3 



(i^^l3,//l3) 



where 



£•3^2 = ^31-82^3 + &32S2dl + 633-82^2 + ^34^2-03 + 6355161 +636i?l 63 

+637B361 + bssBsbs 

C2dz = B31D3C2 + B32D1C2 + B33D2C2 + ^34^3^2 + B35C1C1 + B3QC3C1 

+B37C1C3 + B38C3C3 (A.36) 



631 = - — ^-(Q2)6l4-(0l)624 
2^3 A3 

, n y7 Xi xi Y5 XiYt , n n , p n ^ 

O32 — —J2 -7^ -\ -T^^13-\ ^^(,i>ll + -'>'l3£'16 + -C/03i'15J 

X3 A3 2:3 A3 2:3 A3 

-(Q2)&12 - (gi)^'22 

&33 = -/3 — ^ + —^Du + —^{DuBie + E02B15) - (Q2)6i3 - (gi)6: 

2^3 -'^3 2:3 A3 X3 A3 

2:1 Ai Xt Yv , , 2:1 yt; 

^'34 = —^ + —^{Bi4 + DnBie) + —^Di7 
X3 A3 2:3 A3 2:3 A3 

635 = —^Di2 + — ^Pl2Bl6 + E01B15) - (Q2)6l5 - (Ql)625 
2:3 A3 2:3 A3 

he = -r^ - (Q2)&i6 - (Qi)&26 

2^3 A3 

^37 = -^i?ll + -^i?lli?16-(gi)&27 
2^3 A3 2:3 A3 

b28 = --|^-(Qi)628 (A.37) 

2^3 A3 



'23 
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Here we have used the notation 

0i=F3 + ^f, Q. = F, + f,^f + ^f (A.38) 

X3 A3 X3 A3 X3 A3 

The remaining commutation relations do not participate effectively for the algebraic BA. Neverthe- 
less, we will write them below, but without their coefficients which are very cumbersome. 

(F^25,//25)and(F27,/27) 

L»2&3 = X2lb3D2 + X22blD2 + X23Bidi + X2iBid2 + XasSida 

+Ar26-B3'^l + X27B3d2 + X28-B3'^3 + -^^29-6201 + X210B2C3 
-\-x2ll02cl + X21202C3 

(73^2 = a:2id2C3 + a:;2id2Ci + a;23-DiCi + a;24-D2Ci + X2^D3Ci 

+X2f>DiC3 + X27D2C3 + X2SD3C3 + X29B1C2 + X2l(3B3C2 

+x2iibiC2 + x2i2b3C2 (A.39) 

and 

C162 = Ynb2Ci+Yi2b2C3 + Yi3B2Ci+Yi4B2C3 + Yi5biDi 

+YiebiD2 + YnbsDi + Yi8&3£'2 + YigB^di + YiioBidz 
+YinBid3 + Yii2B3di + ^11353^2 + ^11453^3 

C2&1 = yii&iC'2 + ^1263(^2 + V13B1C2 + yiAB3C2 + yi^diCi 

+yi6d2Ci + yndiC3 + yisd2C3 + y^DiCi + yiwD2Ci 

+2/111-D3C1 + yii2DiC3 + yii3D2C3 + yii4-D3C3 (A. 40) 



• 



(i^J^24,//24)and(F37,/37) 



C262 = C3iD2d3 + C32B2C2 + C33D3d3 + £341)2^2 + C35-B3C1 

+ 036^2-03 + C3jDid2 + C3sdiD3 + C3gdiD2 + C3io&3C'i 

+ C311B1C1 + C3i2D3di + C3i3d3Di + 0314^2 -D2 + C3i^D3d2 

+C3wDid3 + C3i762C'2 + caisdiDi + C3iQDidi + C32o-D2<ii 

+ 0321^3^3 + C322biCi + C323(i2-Dl + C324fclC'3 + C325-B3C3 (A. 41) 
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and 



C363 = C4idiD2 + C42Did2 + C43D2di + C44B1C1 + C45d2D3 

+C4663C1 + C4762C2 + C4sD3di + C4gd2Di + C4iQD^d2 

+C41161C3 + C4i2d2D2 + C41361C1 + C4i4Did3 + Cn^diD^ 
+ 0416-02^2 + C417rfl-Dl + CiiaDidi + C4i9i32C2 + C420(i3-Dl 
+ C421^3C3 + C422^3C3 + C423^3Cl + 2424-03^3 + C425^2rf3 (A. 42) 



• 



(FFi2,//l2)and(F^45,//45) 



and 



£'363 = Xaib^Da + X32biD3 + XssBidi + X34Birf2 + ^35^1^3 

+-'^36-83^1 + -^^37-83^^2 + X^sB^d^ + X39S2C1 + X310-B2C3 

+-''^311&2C'l + -^312^2C'3 

C3d3 = a:;3i(i3C3 + a:;32d3C'i + a;33Dici + X34D2C1 + X35D3C1 

+X3eDiC3 + X3tD2C3 + a;38i:'3C3 + X3gBiC2 + X3i(3B3C2 

+X3iibiC2 + a;3i263C2 (A.43) 



C362 = Y2ib2C3 + Y22b2Ci + Y23B2Ci+Y2iB2C3 + Y2rMDi+Y2QbiD2 

+ Y27biD3 + Y2sb3Di + Y29b3D2 + Y21Qb3D3 + Y21lBidi + Y212Bid2 
+ Y2l3Bid3 + Y2uB3di + Y2l5B3d2 + Y2wB3d3 
C263 = y2lb3C2 + y22blC2 + y23BlC2 + y24B3C2 + y25dlCi + y26d2C\ 

+y27d3Ci + y28diC3 + y29d2C3 + ^210^3*^3 + y21lDiCi + y212D2Ci 

+y2i3D3Ci + y2iADiC3 + J/215-D2C3 + J/216-D3C3 (A. 44) 
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